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Abstract. We investigate the nature of the BL Lac candidates in the northern part of the 2BL sample ILondish et al. (2002)^ 

searching for optical variability by means of Rc band photometry with the Asiago 1.8 and Loiano 1.5 meter telescopes during 
years 2002, 2003 and 2004. We also made aperture photometry of the candidates on the plates available on the web from the 
Palomar and UK Schmidt telescopes. No significant variability was detected for the majority of the objects without detected 
radio and/or X-ray emissions, while clear variability was found for the few sources of the sample with strong radio and/or 
X-ray emission, which is a constant feature of the "classical" BL Lacs. Some considerations on the nature of the 2BL sources 
are made, taking into account their most recent data available in proper motions catalogues (SuperCOSMOS, USNO-Bl) and 
in the Sloan Digital Sky Survey. 
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1. Introduction 

, The history of the BL Lacertae objects started in 1968 
' dSch mitt 1968 1 with the discovery of the radio emission 
, from the prototype of the class, the variable "star" BL Lac. 

A few years later, after that some other similar objects 
: have been found (OJ 287, W Com, ON 325 and AP Lib), 
' IStrittmatter et al. 19721 made the first definition of the class. 
, Subsequent observations with X-ray satellites showed that BL 
Lac objects were also remarkable X-ray sources. The first at- 
, tempts to derive a sample of BL Lacs starting from catalogues 
■ of optically variabile sources were unfruitful, while those based 
on Flat Spectrum Radio Sources provided soon a good deal of 
objects, indicating that the radio emission was a better discrim- 
inant. Several samples of BL Lac objects have been collected 
up to now, using different selection criteria (e.g. 1 Jy sample, 
IStickel et al. 1991i Extended Medium Sensitivity Survey, 
IMorrisetaL 19911 IWolter et al. 19941 Rosat-Greenbank 
survey, 'Laurent-Muehleise n et al. 19991 Sedentary Survey, 
IGiommi et al. 1999. DXRBS, ILandt et aL 200U HRX-BL 
Lac survey, IBeckmann et al. 2003> . In all these samples, the 
presence of a non thermal radio or/and X-ray emission was 
required. 

Nowadays, BL Lac objects are considered a subsample of 
the population of the Active Galactic Nuclei, and are gener- 
ally defined as sources with featureless optical spectra, strong 
and fast variability in the radio, optical and X-ray bands, strong 
and variable radio and X-ray emission and often with a high 
degree of polarization. The precise definition of the BL Lac 
class is still a matter of debate: the optical spectrum was de- 
Send offprint requests to: R. Nesci 



fined as "featureless" if the equivalent widths of the emission 
lines were less than 5 A and/or the Call break was less than 
0.25 ( IStocke et al. 199lV This definition was later relaxed by 
[Marcha et al. (1993)| who included in the BL Lac class also ob- 
jects with stronger features (e.g. Call break < 0.40). A recent 
discussion of these criteria was made bv .Landt et al. 20041 who 
pointed out that the relative intensity of the AGN and of the host 
galaxy is a key factor for the optical line equivalent width. 

Recently [Londish et al. (20 02) defined the first "optically 
selected" sample of BL Lac candidates, the 2BL sample, as 
a by-product of the 2QZ sample of QSO (Croom et al. 200ll. 
Two strips of sky were surveyed by the 2QZ project: one cen- 
tered on the celestial equator and one at declination -30. The 
selection of the QSO candidates was made on UK Schmidt 
photographic plates scanned with the APM, searching for 
point-like sources with photographic (u'-b) and (b-r) Quasar- 
like colours. Spectra for these candidates were then obtained 
with the 2dF fiber spectrograph mounted on the 3.9 m AAT 
telescope. BL Lac candidates were selected from these spectra, 
looking for objects in the magnitude range (19.97> b >18.25) 
and appearing featureless. High proper motion sources were 
then rejected as probable DC white dwarfs. 

The remarkable feature of this sample is that the large ma- 
jority of the sources (47 out of 56) have no radio emission 
at the NVSS JCondonetaL 1998> sensitivity level (2.5 mJy) 
and that only 5 sources have been detected on the ROSAT All 
Sky Survey (RASS; Voges et al. 1999, Voges et al. 2000) , all of 
them being radio loud. Only three of these RASS sources were 
already reported in the literature as BL Lacs. The lack of X-ray 
and Radio emission in the large majority of the 2BL sources is 
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puzzling, because it implies an emission model somehow dif- 
ferent from that generally accepted for BL Lacs. 

The Spectral Energy Distribution (SED) of BL Lacs is 
indeed well represented, in the Log(v)-Log(vFv) plane, by a 
double-bell shape. The first bump may be well described by 
a log-parabohc law (ILandau et al. 19861 IMassaro et al. 2004al 
and can be explained by a synchrotron radiation from rela- 
tivistic electrons moving in a jet aligned along our line of 
sight; the second one, from X to y-ray energies, is presently 
less well constrained in shape by the observations and is 
generally attributed to Inverse Compton scattering of low- 
energy photons on the same relativistic electrons. According 
to the peak frequency of the first bump, BL Lac objects 
are further subdivided into HBL (high-energy peaked) and 
LBL (low-energy peaked), following the scheme introduced by 
IGiommi & Padovani (1994)1 

The faintness (actually upper limits) of the radio/optical 
and X-ray/optical ratios of the 2BL objects means that the low- 
energy bump of the SED is much narrower than that of the 
"typical" BL Lacs. There is however no doubt about the ex- 
tragalactic nature of a number of 2BL sources: some of them 
have indeed a definite redshift from faint emission or absorp- 
tion lines. The most notable case is 2QZ J215454. 3-305654 
JLondish et al. 20()4l l. a radio- and X-quiet source which ap- 
peared featureless in the 2dF spectrum but showed a redshift 
z - 0.494 in a high S/N spectrum obtained with the ESO-VLT. 
A thermal origin for the bulk of the emission (e.g. from an ac- 
cretion disk) rather than from a synchrotron process might be 
an explanation, as suggested bv iLondish et al. 20041 this raises 
however the question whether such sources can still be properly 
classified as BL Lacs. 

To further explore the characteristics of the 2BL candi- 
dates we performed a search for optical variability, which 
is one of the defining features of the class. In the optical 
range the LBL sources show a large variability (3-4 mag- 
nitudes), while HBL sources show a moderate variability 
range (about half magnitude, se e e.g. [He idt & Wagner 1996 
IKurtanidze & NikolashviU 20021 IMarchili et aL 200 4 1. It is, 
however, unlikely that the 2BL candidates are of the HBL type: 
actually their plot in the radio-optical and optical-X spectral in- 
dex plane (a^o - gpxJ Londish et al. (2002)| fig. 13) shows that 
all the X-ray undetected sources have aox ^ 1 0, while HBL 
sources have aox ^ 0.9 (peckmann et al. 2003). We expect 
therefore to detect a variation of the order of at least several 
tenths of magnitude if the 2BL candidates behave like "nor- 
mal" LBLs. 

In this paper we report the results of our optical observa- 
tions, covering a time interval of three years (2002-2004), of 
the 2BL sources well observable from northern mid-latitudes 
(the equatorial 2BL sample). Furthermore we retrieved from 
on-line archives (ESO, STScI) a number of the Palomar and 
UK-ESO Schmidt survey plates and derived photographic red 
magnitudes to explore the historical behaviour of these sources. 
Finally, we considered the Sloan Digital Sky Survey Data 
Release 2 (SDSS-DR2, http://www.sdss. org/dr2/l photometric 
and proper motion data, recently availables for all the sources 
of the equatorial sample. 



Table 1. Observation log with the Asiago and Loiano tele- 
scopes. 



Year 


Month, Day 


Observatory 


2002 


February 11,18 


Asiago 


2002 


March 7-10 


Loiano 


2002 


March 19 


Asiago 


2002 


May 13 


Asiago 


2003 


February 24 


Asiago 


2003 


March 22-25 


Loiano 


2003 


April 23 


Asiago 


2004 


February 14-15 


Asiago 


2004 


February 20-24 


Loiano 


2004 


March 23-24 


Loiano 


2004 


March 29-30 


Asiago 


2004 


April 26-29 


Asiago 



In October 2003 Londish revised the 2BL sample in her 
PhD thesis ( [Londish (2003)| l using, for a number of sources, 
updated proper motions from the SuperCOSMOS sky survey, 
5 band photometry from the Sloan Digital Sky Survey Early 
Data Release ( Sto ughton et al. 2002^ , her own JHK and opti- 
cal photometry. These data allowed to recognize a number of 
candidates in the former ( |Londish et al. (2002 ) i list as probable 
white dwarfs (from proper motion, variability or black-body 
spectral energy distribution). We present anyway in this paper 
our observational results for all the sources monitored by us, 
even if some of them are most likely stellar. 

2. Observations 

2.1. CCD observations 

The 2BL sample ( P^ondish et al. (2002)^ consists of 56 ob- 
jects. We found that two sources in the equatorial strip, 2QZ 
J105534.36-012617 and 2QZ J121834.8-011955, both de- 
tected in the NVSS, were already classified as BL Lac objects. 
The former was identified because detected in the RASS with a 
flux of 3.3 X 10 '^er§ i 'cm^^ (Appenzeller et al. 1998^ while 
the latter was known as the optical counterpart of a radio source 
in the PKS catalogue (Condo n etal. 197 61 

Our follow-up observations were made during the years 
2002, 2003 and 2004 in the Rc band with the Asiago (1.8 
m) and Loiano (1.5 m) telescopes, respectively equipped with 
the AFOSC and BFOSC instruments, covering a field of view 
slightly larger than 8 arcmin. Some Asiago observations were 
made in service mode. Weather conditions did not permit a uni- 
form time sampling for all the objects. Table^reports the dates 
of the observing runs. 

Twentynine out the 3 1 sources in the equatorial strip were 
observed. The typical integration time for the esposures was 
10 minutes; two exposures were made for the fainter sources 
and then averaged. The three "radio detected" candidates 
of the equatorial strip (i.e. 2QZ J105534.36-012617, 2QZ 
J121834.88-011956 and 2QZ J142526.20-011826) were ini- 
tially not included in the monitoring list so our sampling for 
them is poorer. Two sources, 2QZ Jl 14221.4-014812 (red- 
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shift 1.276) and 2Q Z J 1 14327.3-005050 (redshift z = 1 .59 1 ; 
ILondish et al. (2002)> , were never observed. 

Differential photometry of the BL Lac candidates was per- 
formed using a set of about ten stars in the field of each source, 
spanning a range of two or three magnitudes around the object 
of interest. These were selected so that 5 stars were nearly as 
faint as the target, to estimate the actual photometric accuracy 
achievable (check stars), while the brighter ones were used to 
put on a common scale the instrumental magnitudes from dif- 
ferent runs (reference stars). The zero point of the magnitudes 
for each field has been established adopting the GSC2 mag- 
nitude for the brightest reference star. A table of the adopted 
reference stars and magnitudes is given in Appendix (Table 
Al). For each frame a photometric radius equal to the aver- 
age FWHM of the reference stars was used. Aperture photom- 
etry was carried out using the DAOPHOT routine of the IRAF 
package. All the 2BL sources showed a point-like PSF in our 
images. 

During the consistency check between the instrumental 
magnitudes obtained with the two telescopes, an appreciable 
systematic difference was found for a few stars: inspection of 
their colours using the public SDSS magnitudes showed that 
this difference was present only for "red" star (SDSS r - i > 
0.7). These stars were then discarded for the magnitude zero 
point definition in the final data reduction. In no field did we 
find more than one "red" star. The typical average accuracy of 
our magnitudes as a function of a star luminosity can be derived 
from Fig.n where the rms deviation of the instrumental mag- 
nitudes for each check star of all our fields is plotted against 
its average magnitude, separately for the Asiago (crosses) and 
Loiano (squares) data sets. The general trend of this plot is of a 
larger dispersion for fainter stars, as expected. Some stars how- 
ever show a very large dispersion, so we considered them as 
variables and discarded them as check stars in the data analy- 
sis. 

To discriminate between true variability of a BL Lac can- 
didate and photometric uncertainty we binned the magnitude 
values of the check stars at steps of 0.4 mag and computed 
for each bin the average rms value and the dispersion of the 
data around this average. Then we assumed as upper limit of 
the "statistical" deviations for each bin its average rms plus the 
dispersion. A smooth line approximating these upper limits is 
shown in Fig. ^ and was adopted as the border-hne for detect- 
ing variability. 

Table |2] gives the source name, the photographic red mag- 
nitude derived from the b and {b - r) values published by 
ILondish et al. (2002)| the average Rq magnitude, rms deviation 
and number of observations for the Loiano and Asiago datasets 
and the border-line value of the photometric uncertainty. 

A first consideration from the data of Table |2] is that, for 
the large majority of the sources, the average magnitudes in 
the Asiago and Loiano datasets are very similar, but the Loiano 
values are systematically slightly brighter, on average by 0.04 
mag. Given that the Loiano and Asiago observations were 
made each year typically at an interval of about one month (see 
Table [2 it is unlikely that we observed our sources always in 
a brighter state at one telescope, so we attribute this difference 
mainly to the presence of a colour equation in the instrumen- 




18 19 20 

R 



Fig. 1. The rms deviation of the faint reference stars as a func- 
tion of instrumental magnitude for the Loiano (squares) and 
Asiago (crosses) data. The continuous line marks the adopted 
boundary between true variability and photometric uncertainty 
(see text). 

tal magnitudes rather than to intrinsic source variability. We 
applied therefore a systematic correction of -0.04 mag to the 
Asiago (the less numerous) data quoted in Table|2lbefore merg- 
ing the two datasets. 

A second consideration is that generally the rms deviations 
of both datasets are smaller than the borderUne derived from 
the comparison stars, highly suggestive of a limited, if any, op- 
tical variability, at variance with the expectations for BL Lac 
objects. 

The resulting light curves are shown in Figs |5] to |6l where 
we plot for each object our magnitudes and error estimates. The 
average value is shown as a dashed line, while the range (+lcr) 
of the distribution of the magnitudes of each source is shown 
by the dotted lines. 

We found large variability only for two sources: 1218-01, 
a radio loud (not X-ray detected) source, and 1 106+00, which 
has a hot (32000 K) thermal spectrum according to the SDSS 
photometry jLondish (2003)t . The latter source is very faint 
(R ^ 20.3, by far the faintest of our sample), and our photome- 
try is rather inaccurate at this flux level (expected rms - 0.14) 
so we are unsure of its variability. 

A comparison of our data with the photographic 2QZ r val- 
ues in Londish et al. (2002)| (column 3 of our Table |2j is not 
straightforward because of the possible presence of systematic 
effects in the magnitude scale for each field, and to the lim- 
ited accuracy of photographic plates at these magnitudes. Both 
APM and GSC2 magnitudes have a formal error of 0.3 mag, 
which may be even larger at r ^ 19.5. We have checked for our 
"bright" (R < 18.5) reference stars the consistency between 
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Table 2. Observed mean Rc magnitudes and rms deviations for our Loiano and Asiago observations. The last column report the 
fiducial border line value, discussed in the text, for each source. 



Source name 




Loiano 


rms 


IN 


Asiago 


rms 


iN 


border-line 


J lUU^j UUi /zo 


1 n m 


1 n r\i 


U.U4 


4 


1 n n^ 


U.Uj 


o 



n n^ 
U.Uj 


J iUzoi J.J— UUUoJU 


1 n 1 T 
ly.l / 


iy.U4 


n no 
U.Uz 


o 


1 n c\i 


U.U3 


c 

J 


n nc 
U.Uj 


o/^v T 1 n^Am A 1 n 1 ^A^Q 

ZV^A J iUjOU/.4+Ul JDJO 


1 ^"7 


lo.j 1 


n 

U.Uj 


c 

J 


1 O ^ 1 

lo.jl 


U.Uj 


c 

J 


n r\A 
U.U4 


1/^7 T 1 A/i ^ 1 n T 1 nnoA i ^ 
zv^Z J iU4Di5'. /+UUzDij 


lo.4J 


1 /I Q 


n no 
U.Uz 


c 

J 


1 o ^n 
Io.jU 


U.U2 


A 

4 


n n/i 
U.U4 


o/^v Tin^^^^ 1 nn^^TO 

ZV^A J lUjJJJ. 1— UUjJJO 


1 n no 


1 n 


n m 
U.U / 


/ 


1 n on 


U.U4 


c 

J 


n nA 
U.Uo 


zv^Z J iUjjj4. J— Ulzoi / 


1 

Io.jU 


1 no 


n ni 
U.Ul 


L 


1 O 1 o 

lo. Iz 


U.Uj 


A 

4 


n m 
U.Uj 


o/^v Til riA/i /I ^ 1 nnm i n 
L\^l^ J i lUo44.j+UUU/ i / 




zU. J J 


n o^ 
U.zj 


c 

J 


on "^n 
zU. Jy 


U.Uj 


11 



n 1 /I 
U.14 


zvjZ J 1 1 j41o.4+UUlU41 


1 


1 1 Q 

lo. lo 


n no 
U.UZ 


"7 


1 Q OA 

lO.ZD 


U.Uj 


4 


n m 
U.Uj 


zv^Z J i 1 jyUU.J— UZU14U 


1 TA 
lo. /O 


1 TO 

lo. /Z 


n nA 
U.UO 


c 

J 


lo. /4 


U.Uj 


4 


n n/1 
U.U4 


oov Ti i/inin ^ nmmA 
zv^Z J 1 14U1U.J— UUzyjo 


1 n m 


1 n oo 
ly.oo 


n nn 

u.uy 




1 n OA 


U.U / 


z: 




n 1 
U. Iz 


or^v Ti 1 /1 1 n 1 nmTjn 
zl^Z J i 141 J / . 1— UUz / jU 


1 T 1 

lo. / 1 


1 0^ 

lO.OD 


n nA 
U.Uo 


/ 


1 o n^ 




1 

1 


n n^ 
U.Uj 


zl^Z J 1 14DZ1.D— UZ4/JO 


1 

lo. J J 


1 /1A 

lo.4o 


n ni 
U.Ul 


c 

J 


1 o ^n 
lo.DU 


n n 1 
U.Ul 




n c\A 
U.U4 


ZV^Z> J 1 14j j4.o+Uu1UZj 


1 Q /I A 

IV. 40 


1 Q /I ^ 
IV. 4j 


n n^ 
U.Uj 


A 




1 O /I G 
1V.40 


U.U4 




z 


U.Uo 


ZV^Z J 1 1 jyUy.D— UZ4Dj4 


1 nn 
lo. / / 


1 AO 
lO.DO 


n c\A 
U.U4 


C 

D 


1 O 1 A 

lo. /4 


U.Uj 




n n/1 
U.U4 


zv^Z J IzUUl j.o+UUUDjz 


1 A^ 


1 

lO.JO 


n no 
U.Uz 


c 
D 


1 O Ayl 

lo.04 


U.U/ 


z 


n n/1 
U.U4 


zv^Z J IzUjjo. 1— UU4zlu 


1 /1A 

lo.4o 


1 AO 

lo.oz 


n no 
U.Uz 




1 O AT 
lo.OJ 


n n 1 
U.Ul 


Q 


n n/1 
U.U4 


zv^Z J 1 zUoU 1 . — UU4z 1 y 


1 OA 

lo.oO 


1 nr\ 
lo.VU 


n no 
U.Uz 


4 


1 o nA 


n n 1 
U.Ul 


z 


n n^ 
U.Uj 


zv^Z J lzloj4.o— Ul lyjj 


1 

lo.zj 


1 A AH 


n 1 A 
U. lo 


c 


1 A O^ 
lO.ZJ 


U.Uj 


z 


n ni 
U.Ul 


zv^Z J izzJjo.U— Ul joiy 


1 0^ 

lo.oD 


1 OA 

lo.Vo 


n ni 
U.Ul 


c 
D 


1 n f\A 


U.Uj 


1 


n n^ 
U.Uj 


zv^Z J 1zj4j /.d— UlZyj J 


1 m 
lo.y / 


1 ni 
lo. / J 


n n^ 
U.Uj 


c 
D 


1 Q TA 
lo. /D 


U.U/ 


1 


n n/1 
U.U4 


ZV^Z J iZj4jJ. /— Ul ioZZ 


1 S S 1 

10.01 


1 S 87 
10.0 / 


U.U4 


A 
4 








U.U4 


2QZ J130009.9-022601 


18.67 


18.77 


0.03 


6 


18.78 


0.01 


3 


0.04 


2QZJ131635. 1-002810 


19.53 


19.44 


0.03 


6 








0.08 


2QZJ13281 1.5+000227 


19.79 


19.71 


0.09 


5 


19.77 


0.05 


3 


0.10 


2QZ J140021.0+001955 


19.30 


19.53 


0.07 


8 


19.51 


0.03 


3 


0.08 


2QZ J140207.7-013033 


20.14 


19.73 


0.08 


8 


19.87 


0.05 


4 


0.10 


2QZ J140916.3-000012 


18.47 


18.41 


0.02 


8 


18.45 


0.01 


2 


0.04 


2QZ J141040.2-023020 


19.27 


18.68 


0.01 


6 


18.73 


0.03 


6 


0.04 


2QZ J 1425 26.2-0 11 826 


19.47 


19.24 


0.01 


2 


19.31 


0.01 


2 


0.08 



GSC2 and our CCD magnitudes, finding always agreement 
within the formal accuracy of the GSC2. Consistency between 
our magnitudes and those given in the on-line APM catalogue 
(http://www.ast.cam.ac.Uk/~apmcat/i was good (i.e. within the 
formal error) for stars brighter than R ~ \9, while strong differ- 
ences were found for the fainter sources. We decided therefore 
to consider as real only a difference greater than 0.3 mag be- 
tween our values and those derived from Londish et al. (2002)1 
Using the 2QZ r photographic data, the variability of 
1218-01 is confirmed and two more sources, 1055-01 (also 
with a NVSS radio detection) and 1410-02 (radio quiet), are 
found to be varied. The variability of 1234-01, 1400+00, 
1402-01 and 1425-01 is also possible, while we cannot get 
any firm conclusion regarding 1106+00 because only a lower 
limit (r > 20.9) can be derived from .Londish et al. (2002)| 



2.2. Historical photographic data 

For each object we retrieved all the scanned images avail- 
able on-line coming from the POSSl and POSS2 Surveys and 
the UK-Schmidt plates in the red (E or F emulsion) band. 
Differential aperture photometry of the 2BL objects was made 
with the DAOPHOT routine of the IRAF package, using the 



same reference stars selected for the CCD frames: the photo- 
metric radius was taken equal to the average FWHM of the 
bright reference stars. Conversion from photographic magni- 
tudes to our CCD Rc magnitude scale was made building a cal- 
ibration curve for each image by means of the reference stars: 
the CCD magnitude of the 2BL object was determined by inter- 
polation between the nearest experimental points of this curve 
and it had an uncertainty of about 0.2 magnitudes. Only in a 
few plates was the object found to be apparently lower than our 
faintest reference star: in this cases we extended our CCD se- 
quence to encompass the apparent magnitude of the candidate 
BL Lac on the plate. 

The historical (in our CCD-/?c scale) magnitudes so de- 
rived are reported in Table|3lfor each object and epoch: the first 
column gives the source name shortened to the most significant 
digits of RA and DEC (HHMMsDD); subsequent columns give 
our estimated Rc magnitude for each year For some fields we 
found more than one plate in a given year, and their magnitudes 
are reported in subsequent lines. When the source was too faint 
on the plate to be reliably measured we report a lower limit. 

Given the low accuracy achievable for the magnitudes of 
our sources from the Survey photographic plates, we assume 
that variability is believable only if a variation of at least 0.6 
mag is detected. Small amphtude variability of our objects is 
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Fig. 2. Light curves of 2BL candidates: abscissa is in Julian Date -2450000, ordinate is R mag. The triangle on the left of the 
vertical bar indicates the value from Londish et al. (2002) located at arbitrary JD. Black dots are our CCD data with error bars 
derived from the internal consistency of the faint reference stars in the data set. The dashed line is the average value of the source 
in the data set; the dotted lines indicate the +lcr strip (see text). 



therefore not testable in this way. Only eight sources seem to 
be variable: 1036+01, 1053-00, 1055-01, 1141-00, 1218-01, 
1223-01, 1410-02 and 1425-01. All the three radio loud 
sources are variable and two of them (1055-01 and 1218-01) 
varied most significantly, as expected for classical BL Lac ob- 
jects. 



2.3. SDSS and proper motion data 

At the time of writing this paper (October 2004) all the ob- 
jects in the equatorial strip of the 2BL sample (save the ra- 
dio loud source 1055-01) have been observed photometri- 
cally (and several also spectroscopically) by the Sloan Digital 
Sky Survey (SDSS) and publicly available through the Data 
Release 2 (DR2). This release covers a larger sky area and su- 
persedes the Early Data Release (EDR, Stoughton et al. 2002) 
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Fig. 3. Light curves of 2BL candidates; abscissa is in Julian Date -2450000, ordinate is R mag. Symbols as in Fig. 2. 



used by 'Londish (2003)' The DR2 photometric observations 
include 5 bands {u',g',r',i',z') covering the optical spectral 
range from 3500 to 9100 A, so that they allow us to check 
whether the spectral energy distribution is consistent with a 
black body (as expected for a DC white dwarf) or with a power 
law (Fv - A v") which is typical of a BL Lac source. Typical 
values for a range between -1 and -2 for LBL sources (e.g. 
[Vagnetti et al. 2003 1, but may be flatter for HBL objects. 

We performed this test for our sources, applying a correc- 
tion of -0.04 to the published u' magnitudes and +0.02 to the 
z' magnitudes, as indicated in the relevant SDSS web page. A 
correction for the galactic extinction was also made, following 



'Schlegel et al. 1998' which anyway is always small (E^b-v) ^ 
0.055), given the high galactic latitudes of the sources. To com- 
pute the of the fit we assumed an error of 0.03 mag for each 
photometric band, including any systematic effect. The results 
came out to be largely independent of the assumed redden- 
ing correction and of the systematic corrections stated above: 
in no case did a "good" black-body fit switch into a power- 
law one (or vice/versa) including/excluding any of the above 
corrections. There are anyway a few cases with no good fit 
with either spectral shape, which remained "bad" also tuning 
the above corrections. The values of spectral indices range be- 
tween 0.0 and -1.4, which are reasonable for AGNs. Our fits 
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Fig. 4. Light curves of 2BL candidates; abscissa is in Julian Date -2450000, ordinate is R mag. Symbols as in Fig. 2. 



of the spectral shapes with a black body distribution give tem- 
perature values in fair agreement with the previous results by 
iLondish (2003), save 1106+00 which is very hot (~ 130,000 
K). 

We report in Table 0] several quantities for each source: 
source shortened name (column 1), the SDSS r magnitude (col- 
umn 2), the SDSS spectral classification (Sp. CI., column 3 
: Unclass means that the spectrum is too poor according to 
the SDSS team for a classification, while n.a. indicates a not 
available spectrum). In column 4 and 5 we report our best 
fit respectively for the spectral shape given as a temperature 
(in K) in case of a black body (BB) or as a spectral index 



a for a power law (PL). Column 6 gives the of the fit: 
there are 3 degrees of freedom in the fit, so values greater 
than 9 means that the fit is acceptable only at a confidence 
level less than 2.5%. Columns 7, 8, 9 give the source proper 
motion according, respectively, to the SDSS, the USNO-Bl 
(Monet et al. 2003 1 catalogue and the SuperCOSMOS Science 
Archive (SSA, http://surveys.roe.ac.Uk/ssa/i: these motions are 
indicated in units of their rms deviation for the SDSS and SSA 
while for the USNO-Bl the proper motion probability is re- 
ported as given by the CDS, with 9 meaning more than 90%. 

An approximate conversion of the SDSS r magnitudes into 
our GSC2-based Rc magnitudes may be made assuming Rc = 
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Fig. 5. Light curves of 2BL candidates; abscissa is in Julian Date -2450000, ordinate is R mag. Symbols as in Fig. 2. 



r-0.23 (see[Fukugita et al. 1996'l. Given that the zero points of 
our magnitude scale for our fields may have a systematic offset 
of about 0.2 mag, we considered as variable only those sources 
which have a magnitude difference, with respect to our Loiano 
values, > 0.3 mag. Three objects appear to be variable from 
this comparison, 1141-00, 1218-01 and 1425-01, the last two 
being radio loud; all these 3 objects appear variable also in the 
historical survey plates. 



3. Variability properties 

Flux variability at optical frequencies is one of the character- 
istics of the BL Lac objects. The most monitored sources are 
generally LBL sources which show no definite time scale and 
a range of variability of a few magnitudes. HBL sources are 
generally considered to be less variable in the optical; however 
most of the well studied HBL's like the so-called TeV blazars 
(e.g. Mkn 421, Mkn 501, lES 2344+514, lES 1959+650) are 
relatively nearby sources with a strong host galaxy, so that the 
apparent luminosity variation in the optical is reduced by the 
steady stellar component of the integrated flux. As told in the 
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Fig. 6. Light curves of 2BL candidates: abscissa is in Julian Date -2450000, ordinate is R mag. Symbols as in Fig. 2. 



Introduction, however, our sources are not likely to be HBL, 
given the undetection of their X-ray flux. 

We collect in Table |5l our evaluations of the variability 
of each source (column 1) in our data (column 2), with re- 
spect to the photographic 2QZ data (column 3), to the older 
historical data (column 4) and to the SDSS photometry (col- 
umn 5). In columns 6, 7 we report the radio flux in mJy re- 
spectively at the frequencies of 1.4 and 8.4 GHz according to 
ILondish et al. (2002)| and |Londish (2003)| The flag X in the 
last column marks objects excluded from the 2BL sample by 
iLondish (2003), 



The results of our variability search may be summarized as 
follows: 

i) only one source (1218-01, radio-loud) showed large 
variability in our CCD monitoring (we will discuss on it in 
more detail at the end of the section); 

ii) three sources (two of them radio-loud) showed variabil- 
ity when the 2QZ photographic data are compared with our 
CCD data; 

iii) eight sources (including all the 3 radio-loud) showed 
variability when the public digitized Sky Survey plates (span- 
ning a 50 years interval) are measured using our reference CCD 
sequences. 
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Table 3. Photographic magnitudes. Plates up to 1956 are from the POSSl survey with 103aE emulsion, later plates are IIIaF 
from the P0SS2 or UK Schmidt surveys. 



name 1952 1955 1956 1985 1986 1987 1988 1989 1991 1994 1995 1996 1997 1998 



1002-00 


19.0 








19.0 








18.9 




1026-00 


19.6 








19.3 




19.1 








1036+01 


19.9 












18.7 
18.7 








1045+00 




18.2 










18.6 








1053-00 




19.0 






19.4 




19.0 




19.8 




1055-01 




18.8 






18.5 








17.6 




1106+00 




>20.0 






20.1 








20.5 




1134+00 


18.2 










18.1 
18.1 








18.4 


1139-02 


19.1 


















19.2 


1140-00 


>19.8 


















>19.8 


1141-00 


>19.5 


















18.7 


1145-02 


18.5 


















18.5 


1145+00 


19.1 


















19.2 


1159-02 




19.0 






18.8 






18.8 
18.8 






1200+00 




19.0 






18.8 






18.8 






1205-00 




18.5 






18.5 






18.5 






1208-00 




18.5 






18.8 






18.8 






1218-01 




17.6 
18.1 








18.6 






18.5 




1223-01 




18.6 








19.2 






19.2 




1234-01 




18.5 








18.6 


18.9 
19.1 




18.6 




1254-01 






18.9 






18.9 


19.0 
19.0 








1300-02 






19.2 


18.8 




18.9 










1316-00 






>19.8 










19.9 


20.1 




1328+00 






19.5 










19.6 


19.5 




1400+00 




19.9 










19.5 


19.5 






1402-01 




>19.8 










20.0 


20.0 






1409-00 




18.7 










18.6 


18.5 






1410-02 




18.5 










18.8 
19.1 
19.0 


18.9 18.7 
18.6 






1425-01 




20.0 










19.7 


19.4 







iv) three sources showed variability comparing our CCD 
data with those of the SDSS, two of them being radio-loud (the 
third radio-loud source is still not observed by the SDSS). 

This gives a total of eight variable sources out of the 29 ob- 
served, as indicated in Table|5] Notably, none of the 3 sources 
with a spectroscopically confirmed stellar spectrum showed 
variability in our study. 

Of these 8 sources: 

a) three are radio loud (1055-01, 1218-01, 1425-01); the 
first two were also detected in the J, H and K bands by the 
2MASS survey, with magnitudes 16.7, 16.0, 15.3 (1055-01) 
and 15.2, 14.4, 13.5 (1218-01); the last two have also SDSS 
5-band photometry which give unambiguous power law SEDs. 

b) one is a radio-quiet QSO (1036+01) with a power law 
spectrum; 

c) one has a power law spectrum (1141-00) and a faint 
radio detection, but strong proper motion in the SDSS/USNO- 
Bl; 



d) three have a possible hot black body spectrum (1053-00, 
1223-01 and 1410-02) and also a possible proper motion de- 
tection: our power-law fit gives for all these sources val- 
ues much worse then the black-body. 1223-01 was excluded 
from the sample by'Lon dish (2003)] due to its proper motion; 
1410-02 has a strong proper motion in the SDSS/USNO-Bl 
catalogues. 

From our combined variability and spectral shape study, 
therefore, it seems that all the 3 radio loud sources among the 
2BL candidates in the equatorial strip have a variability be- 
haviour similar to "classical" BL Lac objects. We remark that 
only 4 of the sources found to be variable are definitely of ex- 
tragalactic nature. 

What about the 21 non-variable sources (at least at our sen- 
sitivity level)? Twelve of them have already been excluded by 
the 2BL sample by |Londish (2003)1 our reanalysis of the spec- 
tral energy distribution and proper motions confirm that they 
are most likely stellar sources (see Table |5}. Of the remain- 
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Source 


r 


Sp. CI. 


BB Temp. 


PLff 


x' 


proper motion significance 




mag 




K 






SDSS 


USNO-Bl 


SSA 


1002-00 


19.17 


n.a. 


10800 




4.0 


9.81 


9 


12.47 


1026-00 


19.28 


n.a. 


11900 




2.5 


6.08 


9 


5.69 


1036+01 


18.80 


QSO z=1.86 




-0.28 


3.1 


0.50 





0.12 


1045+00 


18.62 


STAR 


11900 




0.8 


6.78 


9 


7.61 


1053-00 


19.42 


n.a. 


10200 




11.2 


2.26 


8 


2.18 


1106+00 


20.46 


n.a. 


130000 




238.3 


0.61 





0.27 


1134+00 


18.43 


QSO z=1.48 




-0.71 


11.9 


0.41 





0.25 


1139-02 


19.18 


GAL? 




-1.42 


5.4 


1.20 





0.40 


1140-00 


20.02 


STAR 


13100 




16.6 


1.95 





3.60 


1141-00 


18.64 


Unclass 




-0.61 


16.4 


21.45 


9 


0.67 


1142-01 


18.20 


n.a. 


7300 




24.8 


0.50 





0.00 


1145-02 


18.67 


Unclass 


11300 




2.3 


1.40 


9 


1.97 


1145+00 


19.57 


n.a. 


11300 




2.7 


1.21 





2.29 


1159-02 


19.03 


Unclass 




-0.39 


16.1 


0.56 


9 


0.50 


1200+00 


19.07 


QSO z=1.65 


6300 




91.0 


0.46 





0.56 


1205-00 


18.74 


n.a. 


7300 




7.1 


11.71 


9 


13.23 


1208-00 


19.08 


n.a. 


12400 




4.5 


2.26 


9 


2.58 


1218-01 


17.56 


Unclass 




-1.33 


0.1 


0.68 


9 


0.58 


1223-01 


19.12 


n.a. 


13100 




20.5 


2.34 


9 


2.55 


1234-01 


19.17 


n.a. 


7600 




30.8 


0.68 





1.45 


1254-01 


19.16 


n.a. 


7600 




6.3 


1.88 





2.42 


1300-02 


18.93 


Unclass 


8000 




11.5 


9.09 


9 


10.48 


1316-00 


19.60 


n.a. 




-0.58 


26.2 


3.11 


9 


9.06 


1328+00 


19.77 


n.a. 


10800 




4.0 


2.44 


9 


2.52 


1400+00 


19.65 


n.a. 




-0.03 


4.4 


0.19 





3.86 


1402-01 


19.88 


STAR 


15100 




45.1 


1.08 


9 


0.43 


1409-00 


18.66 


Unclass 


8800 




4.6 


3.52 


8 


3.84 


1410-02 


18.97 


Unclass 


12400 




16.4 


10.69 


9 


1.13 


1425-01 


19.08 


Unclass 




-0.98 


5.6 


0.11 





0.97 



ing 9, one is a spectroscopically confirmed star (1402-01); 2 
are QSOs (1 134+00 and 1200+00); 2 have a black-body spec- 
tral shape and probable detected proper motion (1 145+00 and 
1254-01); only two have a power-law spectrum (1139-02 and 
1159-02), the second one with a faint radio detection at 8.4 
GHz ( |Londish (2003)t . Finally, 1234-01 has an intervening 
absorber at z=1.06, according to Lon dish (2003)| and is there- 
fore most likely an extragalactic source. 

Only 5 of the non-variable sources are surely or probably 
extragalactic. Summed to the 4 variable ones, we get a total of 
9 AGNs, some of them with unclear nature. We have checked 
that the equivalent width of the emission lines in the sources 
classified as QSO in the SDSS are much larger (~ 40 A) than 
the classical 5 A limit by Stocke et al. (1991) or the extended 
limit of ~30 A by Marcha et al. (1996) and should not be con- 
sidered therefore as BL Lac objects. 

A remarkable case is that of 1218-01. It was identified by 
Condon et al. (1976) as the optical counterpart of the radio 
source PKS B 12 16-010 with a featureless spectrum (Wilkes 
et al. 1984). UBV photometry (Adam 1985) indicates that it 
was even brighter than in our observations, and considering the 
faint state reported by Londish et al. (2002) one can estimate a 
variation range of at least 3 magnitudes. A redshift z - 0.415 is 
reported in the literature after Downes et al. (1985) who quote 



a private communication never confirmed by any subsequent 
published line identification. Radio images from VLBA cali- 
brator survey show a compact bright core (Beasley et al. 2002). 
It is therefore a very typical bright BL Lac object included in 
the 2BL sample because of an occasional faint state. 

4. Conclusion 

The main problem we faced on in this work is how to identify 
and classify BL Lac objects. Of course, the mere observation 
of a featureless optical spectrum is not sufficient although very 
specific of this class of AGNs. Detection of a significant vari- 
ability can be a further request, however it could be difficult to 
be noticed because of the possible occurrence of quiescent pe- 
riods and of the randomness of observation sampling. Although 
the failure in detecting a variability does not necessarily imply 
that a source is not a BL Lac objects, the fact that we observed 
significant brightness changes mainly in radio loud sources in- 
dicates that the search for variations can be a useful strategy to 
confirm their nature. 

It is an open problem if a population of "radio quiet", i.e. 
with a radio-optical spectral index less than 0.2 (Giommi & 
Padovani 2004) BL Lac objects does really exist. At present, 
evidence for such a type of AGNs is elusive and the 2BL 
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Table 5. Summary of variability detection and radio flux densities of the monitored sources. Sources with the flag X were 
rejected by Londish (2003) because of high proper motion. 



Source 


Our inonitoring 


Londish 


On line plates 


SDSS 


Radio data Fl^g 












FilAGHzY FiSAGHzf 












mjy mjy 


1002-00 


no 


no 


no 


no 


X 


1 026-00 


no 


no 


no 


no 


X 


1036+01 


no 


no 


yes 


no 


< 0.1 


1045+00 


no 


no 


no 


no 


X 


1W»^»J WW 


no 


no 


yes 


no 


< 0.1 


1055-01 

± W^ ^ W i 


no 


yes 


yes 
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1106+00 


possibly 


7 


no 


no 


< 0.1 


1 1 34+00 


no 


no 


no 


no 




1139_02 


no 


no 


no 


no 




1 140-00 


no 


no 


no 


no 


X 


1141-00 


possibly 


possibly 


yes 


yes 


3.6 


1 145-02 


no 


no 


no 


no 


X 


1 145+00 


no 


no 


no 


no 




1159_02 


no 


no 


no 


no 


0.3 


1200+00 

± ^WW T WW 


no 


no 


no 


no 




1205-00 


no 


no 


no 


no 


X 


1 208-00 

1 j—WD WW 


no 


no 


no 


no 


X 


1218-01 


yes 


yes 


yes 


yes 


244 


1223-01 


no 


no 


yes 


no 


X 


1234-01 


no 


possibly 


no 


no 




1254-01 


no 


no 


no 


no 




1300-02 


no 


no 


no 


no 


X 


1316-00 


no 


no 


no 


no 


X 


1328+00 


no 


no 


no 


no 


X 


1400+00 


no 


no 


possibly 


no 


0.3 X 


1402-01 


possibly 


possibly 


no 


no 


<0.1 


1409-00 


no 


no 


no 


no 


X 


1410-02 


no 


yes 


yes 


no 




1425-01 


no 


possibly 


yes 


yes 
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1. Flux density from NVSS. 



2. Flux density from Londish et al. (2002), Londish (2003) 



sample must be considered partially contaminated by other 
types of objects, mainly of stellar nature. We tried to establish 
what should be the SED of a source to be properly considered 
as a "radio quiet" BL Lac. It is known that the synchrotron 
emission in many LBL objects covers the frequency interval 
from radio to the optical-UV: the energy distribution is gen- 
erally well represented by a log-parabola (.Landau et al. 1986 
IMassaro et al. 2004ab . also combined with a power law on the 
low frequency side (IMassaro et al. 2004b> : 

F(v) = Av"f''+''^°^<'+'''''''^l (1) 

Fig.Qshows two examples of SEDs of two radio loud BL 
Lac objects in the 2BL sample. The peak frequency of the 
brighter (1218-01) is around 10'^ Hz while that of 1055-01 
is in the optical at about 4x10'^ Hz. Of course they are based 
on not simultaneous data and therefore must be considered only 
indicative of the actual distributions. Note that if the power-law 
branch of 1055-01 would be only moderately flatter than the 
one plotted in Fig.0 its radio flux density would be lower than 
5 mJy, making the source only marginally detectable in many 
radio surveys. According to this picture a characteristic of most 
2BL sources is that their synchrotron peak must be at frequen- 



cies higher than lO'** Hz, which correspond to a flat optical 
spectrum. Indeed we see from Tablel^that all the sources with a 
PL spectrum, with the only exception of 1 139-02, have a spec- 
trum flatter than that of 1218-01 and of typical BL Lacs, with 
best fit spectral indices greater than -1.0. Consequently, the 
frequency of their synchrotron peak must be in the optical-UV 
range. We expect therefore that these sources, dilferently from 
1 055 -0 1 and 1218-01, would not be very bright in the near IR: 
in particular, Rq - K colour indices could be smaller than typi- 
cal BL Lac values (i.e. between 3 and 4, e.g. .Nesci et al. 20031 . 
making them too faint to be detectable in the 2MASS survey, 
which has a K-band limit of 14.5 mag (Skrutski eetaL 1995> . 
In any case, near IR follow-up of these sources would be very 
helpful to clarify their nature. 

The detection of sources with featureless optical spectrum, 
without a radio and X-ray counterpart, can be a tool to select 
candidate BL Lacs, however its efficiency cannot be high be- 
cause of contamination from other galactic and extragalactic 
objects. A population of "optically" selected BL Lac objects, 
with characteristic similar to the 2BL sample, should have pe- 
culiar spectral properties implying different emission, geomet- 
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Fig. 7. The Spectral Energy Distributions (SED) of two ra- 
dio loud BL Lac objects in the 2BL sample: 1218-01 (filled 
circles) and 1055-01 (crosses) derived from our photometry, 
2MASS and other Hterature data. Optical points of 1218-01 
are also representative of its variability. Short dashed lines are 
the SEDs modelled with eq.(l) and the long dashed line corre- 
sponds to a flat radio spectrum at a flux density of 5 mJy. 

rical and physical conditions with respect to normal BL Lacs. A 
small number of sources could actually show these characteris- 
tics, and further multiwavelenght follow-up of this sample is re- 
quired to unambiguously pick up them. Our variability search, 
in fact, indicated that a fraction of 2BL sources are not genuine 
AGNs. We expect that, in addition to variability studies like 
that presented in this paper, polarimetric measurements can be 
another useful selection tool as already successfully applied to 
radio sources by Impey & Tapia (1988), who discovered new 
optical faint blazars at high redshift. 
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APPENDIX A 

In this Table we give, for each monitored source, the 
adopted R magnitude of the reference star, used to bring onto 
a common scale the instrumental magnitudes obtained in the 
different runs from the Asiago and Loiano telescopes. 

Column 1 is the shortened source name, column 2 the 
star identification in the GSC2 catalogue, column 3 the Right 
Ascension, column 4 the Declination and column 5 the adopted 
magnitude as reported in the GSC2. 
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Table 6. Standard stars used for the magnitude scale definition 
of each observed field 



Field GSC2name RA(2000) DEC(2000) mag 
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20 
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15. 
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.83 


15. 


.83 
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11 


06 
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15. 


.56 
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11 


34 


24. 
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09 
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16. 


.25 
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11 


38 


55. 
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03 
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16. 
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11 


40 
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01 
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